print('creating spheres')
from calculate_domain_size import *

cx = fdtd_domain["cell_center_x"]
cy = fdtd_domain["cell_center_y"]
cz = fdtd_domain["cell_center_z"]
# print(number_of_spheres)
for ind in np.arange(number_of_spheres):
    from initialize_fdtd_material_grid import material_3d_space
    # distance of the centers of the cells from the center of the sphere
    distance = np.sqrt((spheres[ind]["center_x"] - cx)**2 
            + (spheres[ind]["center_y"] - cy)**2
            + (spheres[ind]["center_z"] - cz)**2)
    I = distance<=spheres[ind]["radius"]
    material_3d_space[I] = spheres[ind]["material_type"]
cx = 0
cy = 0
cz = 0